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Abstract 



Based on the Vlasov-Maxwell equations describing the 
self-consistent nonlinear beam dynamics and collective 
processes, the evolution of an intense sheet beam prop- 
agating through a periodic focusing field has been stud- 
ied. In an earlier paper |Jl|] it has been shown that in 
the case of a beam with uniform phase space density the 
Vlasov-Maxwell equations can be replaced exactly by the 
macroscopic warm fluid-Maxwell equations with a triple 
adiabatic pressure law. In this paper we demonstrate that 
starting from the macroscopic fluid-Maxwell equations a 
nonlinear Schroedinger equation for the slowly varying 
wave amplitude (or a set of coupled nonlinear Schroedinger 
equations for the wave amplitudes in the case of multi-wave 
interactions) can be derived. Properties of the nonlinear 
Schroedinger equation are discussed, together with soliton 
formation in intense particle beams. 

1 INTRODUCTION 

Of particular importance in modern accelerators and 
storage rings operating at high beam currents and charge 
densities are the effects of the intense self-fields produced 
by the beam space charge and current on determining de- 
tailed equilibrium, stability and transport properties. In 
general, a complete description of collective processes in 
intense charged particle beams is provided by the Vlasov- 
Maxwell equations for the self-consistent evolution of the 
beam distribution function and the electromagnetic fields. 
As shown in [jj]] in the case of a sheet beam with con- 
stant phase-space density the Vlasov-Maxwell equations 
are fully equivalent to a warm-fluid model with zero heat 
flow and triple-adiabatic equation-of-state. 

In the present paper we demonstrate that starting from 
the hydrodynamic equations, and using the renormalization 
group (RG) technique [0, || |], |[| a nonlinear Schroedinger 
equation for the slowly varying single-wave amplitude can 
be derived. The renormalized solution for the beam den- 
sity describes the process of formation of periodic holes in 
intense particle beams. 

2 THE HYDRODYNAMIC MODEL 

We begin with the hydrodynamic model derived in [Q] 
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Here g(x;s) — n(x;s)/N and v(x;s) are the normal- 
ized density and the current velocity, respectively, G(s + 
S) = G(s) is the periodic focusing lattice coefficient, 
v 2 -, = 3Po/2^q is the normalized thermal speed-squared, 
and Pq/uq = N 2 /12A 2 is a constant coefficient [|l|], where 
N is the area density of sheet beam particles, and A is the 
constant phase-space density. Moreover, ip(x; s) is the nor- 
malized self-field potential 
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where (j>(x; s) is the electrostatic (space-charge) potential, 
nib and e b are the rest mass and charge of a beam particle, 
and /3b and 7b are the relative particle velocity and Lorentz 
factor, respectively. Finally, the quantity K is the normal- 
ized self-field perveance defined by 
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In what follows the analysis is restricted to the smooth 
focusing approximation 



G(s) = G = const, 
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and assume that there exist nontrivial stationary solutions 
to (2^1) in the interval x € (— x^~\ x^), and that the 
sheet beam density is zero (g = 0) outside of the interval. 
The change of variables 

^ ^ = ip - Gx^x 
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enables us to rewrite ( 2. 1 ) in the form 
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Clearly, the system (2.4) possesses a stationary solution 
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, Vn = 0, *n = - 
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Here, the uniform density go is normalized according to 
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3 RENORMALIZATION GROUP 
REDUCTION OF THE 
HYDRODYNAMIC EQUATIONS 

Following the basic idea of the RG method, we represent 
the solution to equations ( [O] ) in the form of a standard 
perturbation expansion [|[ in a formal small parameter e as 
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Before proceeding with explicit calculations order by order, 
we note that in all orders the perturbation equations acquire 
the general form 

dg n dv n 

-fr + Q0~^ = On, 
-to+2 go v T — = -— + n , (3.2) 



where the functions a n (£; s) and /?„(£; s) involve contri- 
butions from previous orders and are considered known. 
Eliminating v n and $> n , it is possible to obtain a single 
equation for g n alone, i.e., 
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It is evident that in first order a\ = j3\ = 0. Imposing the 
condition 

Veo 

J d£ei(fc*)=0, (3.4) 
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which means that linear perturbation to the uniform sta- 
tionary density go should average to zero and not affect the 
normalization properties on the interval (0, + 
we obtain the first-order solution 

0i(£,l s ) = X! Ane* Xm ^ ;s \ Xmtes) = UJ m S + mat;. 

(3.5) 

Here, A m are constant complex wave amplitudes, and the 
following conventions and notations 

G 

U-m = -U m , CT = 77' A- m = A* m . (3.6) 

K 

have been introduced. Moreover, the discrete mode fre- 
quencies uj m are determined from the dispersion relation 
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In addition, the first-order solution for the current velocity 
and for the self-field potential can be expressed as 
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In obtaining the second-order perturbation equation 
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(3.3), we note that 
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Thus the second-order solution for the density £>2(£; s) is 
found to be 
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Having determined Qi, the second-order current velocity 
U2(C; s) can be found in a straightforward manner. The re- 
sult is 
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In third order, the functions 03 and j3% entering the right- 
hand-side of equation ( 33 ) can be calculated utilizing the 
already determined quantities from the first and second or- 
ders, according to 
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It is important to note that the right-hand-side of equation 
(3.3) for g3 contains terms which yield oscillating terms 
with constant amplitudes to the solution for 03. Apart from 
these, there is a resonant term (proportional to e lXm ^' s ') 
leading to a secular contribution. To complete the renor- 
malization group reduction of the hydrodynamic equations, 
we select this particular resonant third-order term on the 



right-hand-side of equation (3.3). The latter can be written 

as 
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Some straightforward algebra yields the solution for 
£>3(£; s) to equation ( |3.3| ) in the form 



L>:>, 



(3.20) 



m/0 



where the dots stand for non-secular oscillating terms. 
Moreover, the amplitude 7- > m (^;s) is secular and satisfies 
the equation 
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where the operator C m is defined by 
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We can now construct the perturbative solution for g up 
to third order in the small parameter e. Confining attention 
to the constant stationary density go and the fundamental 
modes (first harmonic in the phase Xm), we obtain 

0(6 s) = 9o + eJ2 [An + e 2 V m {^ s)]e l ^ s \ 

(3.23) 

Following the basic philosophy of the RG method, we in- 
troduce the intermediate coordinate X and "time" S and 
transform equation ( 3.23[ ) to 
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Note that the transition from equation ( 3.23 ) to equation 
( 3.24 ) can always be performed by enforcing the constant 
amplitude A m to be dependent on X and S, which is in 
fact the procedure for renormalizing the standard perturba- 
tion result. Since the general solution for s) should not 
depend on X and S, by applying the oper ator C m (X; S) 
[which is the same as that in equation (3.22) but with 
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where we have dropped the formal parameter e on the right- 
hand-side. Since the above equation should hold true for 
any choice of X and S, we can set X — £ and S = s. 
Thus, we obtain the so-called proto RG equation | 
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Introducing the new variable 
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and neglecting the second order derivatives d 2 /ds 2 and 
d 2 /dsd( m , we finally arrive at the RG equation for the ra- 
th mode amplitude 
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4 THE NONLINEAR SCHROEDINGER 
EQUATION FOR A SINGLE MODE 



Equation (3.28) represents a system of coupled nonlinear 
Schroedinger equations for the mode amplitudes. Neglect- 
ing the contribution from modes with k ^ m, for a single 
mode amplitude A m , we obtain the equation 
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It is easy to verify that T rn is always positive. In nonlinear 



optics equation (4.1) is known to describe the formation 
and evolution of the so-called dark solitons [Q]. In the case 
of charged particle beams these correspond to the forma- 
tion of holes or cavitons in the beam. Since the renormal- 
ized solution for the beam density s) can be expressed 
as 

0(6 «) = Qo + An(6s)e lXm(C;s) . (4.3) 
these holes have periodic structure in space £ and "time" s. 



5 CONCLUDING REMARKS 

Based on the renormalization group method, a system of 
coupled nonlinear Schroedinger equations has been derived 
for the slowly varying amplitudes of interacting beam- 
density waves. Under the approximation of an isolated 



wave neglecting the effect of the rest of the waves, this sys- 
tem reduces to a single nonlinear Schroedinger equation 
with repulsive nonlinearity. The latter describes the for- 
mation and evolution of holes in intense charged particle 
beams. 
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